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Abstract 

We numerically solve the (3+l)-dimensional relativistic hydrodynamical equation 
with a bag-model equation of state, which is one of the simplest models for the first 
order phase transition. Based on the numerical solution, we discuss the space-time 
evolution of the hot fluid produced in ultra relativistic nuclear collisions. Especially, 
the space-time structure of the two-phase co-existing region during the quark-gluon 
plasma-hadron phase transition is analyzed in detail. Finally, we investigate the 
change of the mass spectra of J/tp, ij c and <j) as possible experimental signatures of 
the phase transition. 
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1 Introduction 



The QCD phase transition and quark-gluon plasma is currently one of the 
hottest topics in high energy nuclear physics. [1] The recent high performance 
Lattice simulation suggests that the order of the QCD phase transition in 
the real world, i.e. 2 massless quarks and 1 massive quark case, is first or- 
der. [2] If the QCD phase transition is of the first order, the latent heat will 
be released at the phase transition temperature T c and an isothermal Quark- 
Gluon Plasma(QGP)-hadron co-existing region will appear in the case that the 
change of the whole system is slower than the microscopic relaxation time. If 
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such a scenario will be realized, the space time volume of the T = T c region 
can become significantly large, and physical quantities which are characteristic 
to the phase transition temperature T c will be available as the experimental 
signatures of the phase transition. The purpose of this paper is to investigate 
how large a region with T = T c will be produced at the CERN SPS experi- 
ment and how clearly the hadronic mass shift at finite temperature can work 
as possible signatures of the phase transition. 

As an equation of the state, we here adopt the bag model which is one of 
the simplest models for the first order phase transition. The hydrodynam- 
ical model for the QGP fluid has already been discussed in many papers 
since Bjorken first introduced the simple scaling model in 1+1 dimensional 
expansion [3]. We discussed 1+1 dimensional expansion of the viscous fluid [4] 
and 3+1 dimensional expansion with phase transition [5]. Based on the hy- 
drodynamical model, we also discussed electromagnetic probe of QGP [6]. As 
concerns the treatment of the first order phase transition, Alam et al. and 
Kajantie et al. [7] focused their discussion only on the transverse expansion 
based on the scaling ansatz. Sollfrank et al. [8] discussed a smeared equation 
of state with a small temperature width. In order to focus our discussion to 
the T = T c region, we solved 3+1 dimensional expansion keeping an equation 
of state with the exact first order phase transition. 

Experimentally, at SPS energy, net baryon number is not negligible even in 
the central region. In order to take correctly account of the baryon number 
conservation law, we must solve the coupled equations of hydrodynamics and 
baryon current conservation [9]. However, if we adopt baryon number indepen- 
dent of an equation of state, the dynamics of the four velocity is free from the 
baryon conservation law and the role of the chemical potential is only for the 
calculation of the baryon number at the final stage[10]. Therefore, we may say 
that the baryon free fluid model works well enough in the discussion of space 
time evolution of the fluid with the simple bag model. 

The mass shift of the hadrons in the hot medium is believed to be a promising 
detectable signature of the high temperature region. Hioki et al. proposed J/ip 
and r] c [11] and Asakawa et al. proposed as the possible probe [12]. Following 
their discussions, we apply our hydrodynamical model to the evaluation of the 
mass spectrum of these particles. 



2 Hydrodynamical Equation 

Assuming cylindrical symmetry to the system, we introduce cylindrical coor- 
dinate r, r], r and instead of the usual Cartesian coordinate x, y, z and t, 
where, 
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t — T COsh T], 

z = t sinh 77, 
x = r cos 0, 
y = r sin 0. 



(1) 



The four velocity is given as, 



[/* = JJ T cosh 77 -(- C/ 77 sinh 77, 
£/ z = U T sinh r/ + U 11 cosh 77, 
[/* = [T cos + [^ sin0, 
f^ = [/ r sin0- U* cos0. 



(2) 



Taking into account U^U^ = 1 and cylindrical symmetry, can be repre- 
sented with two variables Y L and Ft, 

U T = cosh Ft cosh(>L — 77), 
U v = cosh F T sinh(F L — 77) , 
[/ r = sinhr T , 

f/^ = 0. (3) 
As is well known, the relativistic hydrodynamical equation is given as, 



where E and P are energy density and pressure, respectively. Usually, thermo- 
dynamical quantities such as E and P are treated as local quantities through 
temperature T = T{x^). 

In order to solve the above hydrodynamical equation, we need an equation of 
state. In this paper, we adopt the bag model equation of state as the simplest 
model for the QCD phase transition of first order. In the bag model, ther- 
modynamical quantities below the phase transition temperature correspond 
to those of a massive hadronic gas. On the other hand, in the high tempera- 
ture phase, thermodynamical quantities are given by those of massless QGP 
gas with bag constant B. In this paper, we assumed pions and kaons in the 
hadronic phase and u-, d-, s-quarks and gluons in the QGP phase. Putting the 
phase transition temperature, T c = 160 MeV, into the condition of pressure 
continuity, we can obtain the bag constant as B = 412 MeV/fm 3 . 



= 0, 



(4) 



and the energy momentum tensor of perfect fluid is given by, 



T» v = EU»U V - P{g^ - U»U V ). 



(5) 
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At the phase transition temperature T c , pressure P is continuous but other 
quantities such as energy density E and entropy density S have discontinuity 
as a function of temperature. We parameterize these quantities during the 
phase transition by using the volume fraction A [11]. We assume that ther- 
modynamical quantities at T = T c are the functions of the space-time point 
through A(x M ) and given as, 



E{\) = XE qgp (T c ) + (1 - A)£ had (T c ), 

S(X) = \S qgp (T c ) + (1 - A)S had (T c ), (6) 

where < A(a+) < 1. This parameterization enables us not only to solve the 
hydrodynamical equation easily but also to estimate hadronic fraction in the 
phase transition region for later use. 



3 Space-time evolution of co-existing region 



Assuming that local equilibrium is achieved at 1 fm later than the collision 
instance, we give initial conditions of the hydrodynamical model on the r = 
t = 1 fm hypersurface. Initial entropy density distribution is given as, 

S(r , r,) = S(T ) exp(- (l ^~^ (M - rj ) - ( ^f^9(r - r )), (7) 

2 ■ a v z 2 ■ a r z 

and parameters we used in this paper are summarized in table 1. 



Table 1. Initial Parameters 





T (MeV) 


Vo 


a v 


r (fm) 


cr r (fm) 


S + Au 


209 


0.7 


0.7 


(32) (V3) _ L0 


1.0 


Pb + Pb 


190 


0.7 


1.0 


(197)(i/3) _ L0 


1.0 



Putting freeze-out temperature as T f = 140 MeV for S+Au and T f = 138 
MeV for Pb+Pb, respectively, we can reproduce hadronic distributions of the 
recent experimental results. Figure 1 and fig. 2 show the pseudo-rapidity dis- 
tribution of charged hadrons and the transverse momentum distribution of 
neutral pions, respectively. Data for S+Au 200 AGeV collisions are obtained 
by WA80[13,14] and data of Pb+Pb 158 AGeV collisions are obtained by 
NA49[15] (rapidity distribution of charged hadrons) and WA98 [16] (trans- 
verse momentum distribution of neutral pions). 

Figure 3 displays the space time evolution of the temperature distribution in 
the Pb+Pb case. The isothermal region at phase transition temperature T = 
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T c appears very clearly. The structure of temporal evolution does not differ 
essentially from the case of smooth phase transition in a small temperature 
region [5,8]. By virtue of the use of the volume fraction A at phase transition 
temperature, we can pick up the volume element with the temperature at just 
the phase transition point. The temperature profile function (fig. 4) shows the 
significant contribution of the space-time volume at T = T c . 



4 Change of the mass spectrum 



As a signal of the hot hadronic region, shifted J/i/j mass-spectrum and i] c 
were proposed by Hioki et al.[ll] and <fi by Asakawa et al.[12] In these papers, 
authors discussed possible signals based on a scaling hydrodynamical model. 
We apply our (3+1) -dimensional hydrodynamical model to these calculation 
and investigate numerically how clear the shifted hadronic mass spectra appear 
in the present CERN SPS case. 



According to ref. [11], the invariant mass distribution of leptonic pair decay is 
given by 

where <£> is the temperature profile function which we have already calculated 
in the previous section. In ref. [11], Hioki et al. evaluated shifted mass M(T) 
and decay width IV(T) based on a c-c potential model and we also use their 
result here. Charmonium density n(x^) was estimated by the kinetic equation 
in ref. [11], and we adopt Maxwell-Boltzmann distribution for n(x^), that is, 

'MT\i , M, 



/mN {MT\* , M. 
" (r) = 3 hd exp(--), 



for simplicity. However, charmonium mass is much larger than the typical 
temperature of the fluid, above n(T) is almost constant and final results were 
dominated by the temperature profile function <3>(T). Figure 5 and 6 show the 
possible J/ip and r) c spectrum for the CERN SPS experiment in comparison 
with the results of ref. [11]. In fig. 5, "critical" and "thermal" correspond to 
the yields of T = T c and Tf < T < T c , respectively. "Cold" in fig. 5 is given 
by, 

1 da Nr lT (T = 0) e/ „ /m\ \ 
-— = — ^ ^{M - M(T)) 

a (i M 1 total 

which corresponds to the naive thermal contribution from the freeze-out hy- 
persurface. Note that, in order to make comparison clear, the results of ref. [11], 
which are the (l+l)-dimensional calculation, were so normalized in fig. 5 and 
fig. 6 as to give the same "cold" contribution to our results. In fig. 5 and 
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fig. 6, quantitative difference between S+Au and Pb+Pb is consistent with 
the difference in the temperature profile functions (fig. 4). 

The yield of J ftp from T c is several orders smaller compared to the cold and 
thermal J/ip. However, the situation of i] c is much better than J/ip: in our 
results, r\ c from the T = T c region is almost the same order as the others. 
Hence, r] c will be a better signature of T = T c than J/ip. These results are 
consistent with ref. [11]. Because of the very small decay width of J/ip and 
7] c , changes of the mass seem very clear and easy to detect. On the other 
hand, their life time is much longer than the hot fire ball and the number of 
particles which decay in the hot medium is very small. Therefore, shifted r] c 
and J /if) are hardly detectable in the present experimental situation, but in 
future experiments with high precision, these can work as good signals. 

Asakawa et al.[12] started their discussion from the number of particles which 
decay in unit time and unit phase space which is given as, 



diV = ^^(T)e-^d 4 xd 3 P , (8) 

where is the degeneracy of and 7 is the Lorentz factor of a moving particle. 
They imported the results of the QCD sum rule for the decaying constant in 
finite temperature, T(T). From (8) we obtain 



d4^ = W // ^u(m(M,m,(T))e^ u ^d 4 - d ^ d ^ d <P> (9) 

where F^M, m^(T)) is the normalized smearing function to take account of 
experimental resolution. F^M, m^(T)) is assumed to have a Gaussian form, 



F (M,m (T)) = — L= e -(^Cr))W (10) 

where a is put as 10 MeV. With use of the above formula, we evaluate the 
spectrum based on our hydrodynamical model. Figure 7 shows the spectrum 
from S+Au and Pb+Pb collisions at CERN SPS energy. The dashed line in 
Fig. 7 stands for the contribution of T = T c and the dotted line for the from 
the hot hadronic region. Quantitative difference between yields in S+Au(a) 
and Pb+Pb (b) are almost proportional to the difference in the temperature 
profile functions(fig. 4). Comparing to S+Au, the ratio of the peak around 959 
MeV to the one around 1020 MeV is large in the Pb+Pb case. 

These curves show clear difference as we expected, but the accumulated result 
(solid line) does not show distinguishable double peaks. In the accumulated re- 
sults, we can find clear double peaks originating from the cold component and 
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hot hadronic contribution, however it would not mean directly the existence 
of the QCD phase transition. 



5 Concluding remarks 

We solved a (3+l)-dimensional hydrodynamical model with first order phase 
transition model(bag model). We can parameterize our hydrodynamical model 
so as to reproduce well the present experimental hadron spectrum. Based on 
the numerical solution we discussed the space-time structure of the isothermal 
region at phase transition T = T c , which appears as a sharp peak in the 
temperature profile function. Applying our numerical solution, we evaluated 
mass spectra of J/ip, i] c and <fi which have been proposed as possible signals 
for the existence of an isothermal T = T c region. According to our numerical 
investigation, if the experimental resolution is as good as a = 10 MeV, we may 
have a chance to find the double-peak of the <p spectrum, but in the present 
CERN SPS Pb+Pb experiment, T = T c peak is indistinguishable from the 
hot hadronic contribution. 

We are indebted to Prof. I. Ohba and Prof. H. Nakazato for their helpful 
comments and their continuous encouragement. We thank many discussions 
between Dr. H. Nakamura, Dr. T. Hirano and other members of Waseda 
Univ. high energy physics group. We also thank Dr. T. Peitzmann and Dr. S. 
Nishimura for discussions about the recent WA98 results. Numerical analyses 
have been done with workstations of Waseda Univ. high energy physics group. 
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Figure Caption 

Figurel:(a) Pseudo-rapidity distribution of charged hadrons in S+Au 200 
AGeV collision in comparison with the data of CERN WA80, and (b) rapidity 
distribution of charged hadrons in Pb+Pb 158 AGeV collision in comparison 
with the data of CERN NA49. 

Figure2: Transverse momentum distribution of neutral pions in S + Au col- 
lision 200 AGeV collision in comparison with the data of CERN WA80(a), 
Pb+Pb 158 AGeV collision in comparison with the data of CERN WA98(b). 

Figure3:The space-time evolution of the temperature-distribution of the Pb+Pb 
158 AGeV collision. Initial temperature is 190 MeV, critical temperature is 
160 MeV and freeze-out temperature is 138 MeV. The flat area corresponds to 
the mixed phase. From these graphs we can see that the QGP phase vanishes 
at t = 2.0 fm and the mixed phase vanishes at r = 8.5 fm and at r = 11.5 fm 
the temperature becomes lower than the freeze-out temperature everywhere 
in the fluid. 

Figure4:The temperature profile function of S+Au 200 AGeV collision(a) 
and Pb+Pb 158 AGeV collision(b). The sharp peak at critical temperature 
corresponds to the huge space-time size of the mixed phase region which is 
caused by the released latent heat. 
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Figure5:The invariant mass distribution of lepton pairs from mass-shifted 
J /if) in the hadronic phase of "critical" , "thermal" and "cold" in S+Au colli- 
sion^) and Pb+Pb collision(b). The dashed line and a square point with a 
dot stand for the results obtained by Hioki et al. In these figures, to make 
comparison clear, the results of Hioki et al. were so normalized as to give the 
same "cold" contribution to our results. 

Figure6:The invariant mass distribution of photon pairs from mass-shifted r\ c 
in the hadronic phase of "critical" , "thermal" and "cold" in S+Au collision(a) 
and Pb+Pb collision(b). The dashed line and a square point with a dot stand 
for the results of Hioki et al. In order to make comparison clear, the results of 
Hioki et al. were so normalized as to give the same "cold" contribution to our 
results. 

Figure7:The invariant mass distribution of lepton pairs from mass-shifted <fi in 
the central rapidity region in S+Au collision(a) and Pb+Pb collision(b). The 
dashed line stands for the contribution from the mixed phase, the dotted line 
for the hadronic phase, the chained line for the contribution from the freeze- 
out hypersurface and the solid line for the total contribution, respectively. 
The peak at 1020 MeV corresponds to the contribution from the freeze-out 
hypersurface. 
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